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We discuss the production of cosmogenic neutrinos on extragalactic infrared photons in a model of 
its cosmological evolution. The relative importance of these infrared photons as a target for proton 
interactions is significant, especially in the case of steep injection spectra of the ultrahigh energy 
cosmic rays. For an ® cosmic ray injection spectrum, for example, the event rate of neutrinos 
of energy above 1 PeV is more than doubled. 
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I. INTRODUCTION 

The assumption that the ultra high energy cosmic rays 
(UHECR) are nuclei (presumed here to be protons) accel¬ 
erated in powerful extragalactic sources provides a nat¬ 
ural connection between these particles and ultra high 
energy neutrinos. This was first realized by Berezinsky 
and Zatsepin P] soon after the introduction of the GZK 
effect 12 • The GZK effect is the modification of the UHE 
proton spectrum from energy losses by photoproduction 
interactions with the 2.7K microwave background radi¬ 
ation (MBR). In the case of isotropic and homogeneous 
distribution of UHE cosmic ray sources, the GZK effect 
leads to a cut-off of the cosmic ray spectrum below 10^° 
eV. The charged mesons generated in these interactions 
initiate a decay chain that results in neutrinos. Since 
the mesons and muons do not lose energy before decay, 
the high energy end of the spectrum of these neutrinos 
follows the injection spectrum of UHECR, while below 
the interaction threshold it is flat i, i. The neutri¬ 
nos which are produced by photomeson producing inter¬ 
actions of UHECR nuclei are sometimes referred to as 
cosmogenic neutrinos. 

Several calculations of of the fluxes of UHE photome¬ 
son neutrinos were published in the 1970s 13) Hi IS S 0 > 
Hill and Schramm (3, 0 used the non-detection of such 
neutrinos to place an upper limit on the cosmological evo¬ 
lution of the sources of UHECR. The problem has been 
revisited several more times [TS Ell [l3 • 

In 2004 Stanev E3l considered interactions of UHECR 
with photons of the extragalactic infrared and optical 
background (IRB), pointing out that this process gener¬ 
ates non-negligible cosmogenic neutrino fluxes. This sug- 
gestion was quickly followed by a confirmation in Ref. [Ij| 
which emphasized the importance of the IRB as interac- 
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tion ta rget . This idea was further developed in Ref. M 

Ref. [I3| gave an estimate of the cosmogenic neutrino 
flux generated in interactions on the IRB, but did not 
account correctly for the cosmological evolution of the 
infrared background. In this paper we perform a cal¬ 
culation using a realistic empirically based model of the 
cosmological evolution of the spectral energy distribution 
of the extragalactic IR-UV background given in Ref. M 
which will be referred to as SMS05. The aim is to es¬ 
timate correctly the role of these extragalactic photons, 
particularly the infrared photons which are by far the 
most numerous, as targets for UHE proton interactions. 

The paper is organized as follows: in Section H we 
discuss the model of the infrared background and its cos¬ 
mological evolution. In Section HI we describe the cal¬ 
culation. Section IV gives the results of the calculation 
and Section V contains the discussion of the results and 
the conclusions from this research. 


II. COSMOLOGICAL EVOLUTION OF THE 
IR-UV BACKGROUND 

It is now well known that galaxies had a brighter past 
owing to the higher rate of star formation which took 
place. Strong evolution is supported by many observa¬ 
tions relating IR luminosity to the much higher star for¬ 
mation rate at z ~ 1 and to the recent determination that 
most Lyman break galaxies at z 1 are also luminous 
infrared galaxies. In addition to the evolution of galaxy 
luminosity, some increase in galaxy number density is ex¬ 
pected owing to the hierarchical clustering predicted by 
cold dark matter models. However, luminosity evolution 
is the dominant effect and it is difficult to separate out a 
component of density evolution. 

In order to calculate intergalactic IR photon fluxes 
and densities and their evolution over time (or redshift), 
SMS05 performed an empirically based calculation the 
SED of the IBR (infrared background radiation) by using 
(1) the luminosity dependent galaxy spectral energy dis¬ 
tributions (SEDs) based on galaxy observations, (2) ob- 
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servationally derived galaxy luminosity distribution func¬ 
tions (LFs) and (3) the latest redshift dependent lumi¬ 
nosity evolution functions, sometimes referred to as Lilly- 
Madau plots. The SMS05 calculation was an improved 
version of the work presented in Refs. Ea, M and M 
The calculation considers two different cosmological 
evolutions, £(z) baseline and fast, of the infrared emission 
of the type 


£{z) 


{1 + z)^ 

(1 + Zflat)™ 
0 


Zflat < 2: < 6 

z > 6 


( 1 ) 


The baseline evolution model is described by m=3.1 and 
2flat = 1-3, while the fast evolution model uses m=4 and 
2flat = 1- The infrared emission at z > Zflat is constant in 
both models. FigureClshows the number density between 



FIG. 2: Number density of the infrared background at dif¬ 
ferent redshifts calculated by SMS05 in the fast evolution 
model. 



FIG. 1: Number density of the IRB at different redshifts as 
calculated by SMS05 Hra . 


III. THE CALCULATION 

The calculation was performed in two stages: (1) cal¬ 
culation of the neutrino yields from interactions with ex- 
tragalactic infrared photons and (2) a subsequent inte¬ 
gration of the yields to obtain the cosmogenic neutrino 
flux from such interactions. This approach gives us the 
flexibility to easily obtain the neutrino flux using dif¬ 
ferent parametrizations of the cosmic ray emissivity, in¬ 
jection spectra and cosmological evolution of the cosmic 
ray sources. This approach, however, suffers from the 
fact that since the yields are calculated only on the IRB, 
they do not account for the fact that high energy protons 
interact mainly with the much more numerous MBR pho¬ 
tons. 


photon energy of 3.16x10“^ and 1 eV in both models. 
The fast evolution model has higher density in the current 
cosmological epoch as well at the IRB maximum epoch, 
which is around z = 2. The increase of the total IRB 
number density increases by a factor of about 4 from z 
= 0 to z = 2 and decreases at higher redshifts. One 
should note, however, that the cosmological evolution of 
the infrared background density is much slower than that 
of MBR since the current IRB density is accumulated 
from the infrared emission of different sources since z = 
6. Figure|21shows the energy spectrum of the fast infrared 
background at redshifts from 0 to 5. One can see both the 
increase of the total photon density as well as the shift of 
the maximum of the emission to higher energy at higher 
redshifts. In terms of photoproduction interactions on 
IRB this means that lower energy cosmic rays will be 
above the photoproduction threshold at higher redshifts. 

Both figures above show the number density of IRB 
rather than the usual presentation of the energy density. 
Since we are using the infrared background as a target 
for cosmic ray interactions this is the relevant quantity. 


A. Calculation of the Neutrino Yields 

The neutrino yields as a function of the proton energy 
Ep, neutrino energy and the redshift, z, were cal¬ 
culated using the IRB spectra at different cosmological 
epochs, i.e., as a function of redshift, that were provided 
by the authors of Ref. E3- Each of the yield calcula¬ 
tions was performed for proper distances corresponding 
to Az=0.2 using an nM=0.3, Oa=0.7 cosmology as 

D{z) = — / [nM(l + 2 )^-I-IIa] (2) 

iiQ Jp . L Z 

The IRB is considered to be constant during each cos¬ 
mological epoch of duration Az=0.2. 

In this way the matrix element corresponding to dt/dz 
dependence was accounted for in the yields. The yields 
were calculated with the code used in Ref. 113 and the 
photoproduction interaction code SOPHIA . All gen¬ 
erated neutrinos are redshifted by the code to the end of 
the Az epoch. The yields were calculated for redshifts 
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0 < z < 5 and for cosmic ray energies above 10^® eV in 
ten logarithmic bins per decade of energy. 



FIG. 3: Neutrino yields for 10^° eV protons interacting with 
both MBR photons and for eV, 10^® eV and 10®° eV 
protons interacting with IRB photons at 2 = 0, given for 
protons traveling a distance of 1 Mpc. 


Fig.El compares the yields for UHE protons travel¬ 
ing a distance of 1 Mpc and interacting with MBR and 
IRB photons. The yield for 10®° eV protons interacting 
with IRB photons is about a factor of 10 lower than that 
for MBR interactions. This difference is much smaller 
than the ratio of the MBR and IRB total densities, and 
demonstrates that 10®° eV protons interact mainly with 
photons in the higher frequency Wien tail of the 2.7K 
MBR spectrum. 

Protons of energy 10^° eV do not interact at 2=0 with 
MBR photons, but they readily interact and produce 
neutrinos by interactions with IRB photons, as do pro¬ 
tons of energy 10^® eV. Even protons of energy 10^^ eV 
occasionally interact with IRB photons, but their con¬ 
tribution is very small and is neglected in this calcula¬ 
tion. Even for £1“® cosmic ray spectra, the smaller 10^° 
and 10^® eV yields are multiplied by the much higher 
flux of cosmic rays with such energies. This is the basis 
of the significant neutrino (and 7 -ray) production from 
UHECR-IRB interactions. 


B. Integration of the Yields 


shift evolution of flat spectrum radio sources, given as an 
analytic approximation in Ref. |2l| . We use the fast evo¬ 
lution model from SMS05 since it is more consistent with 
the new observations of the Spitzer telescope . 

We normalize the cosmic ray energy flux at E„ = 10^° 
eV to EpdNplEp = 2.5x10-1® cm-®s-isr-i ]|l [2|. 
Since the calculation is extended to energies below IQi® 
eV, the code uses the cosmic ray flux at IQi® eV to cal¬ 
culate the injection spectra at lower and higher energy. 
Therefore, the cosmic ray emissivity above IQi® eV de¬ 
pends on the injection spectrum. The injection spectrum 
itself is used as a free parameter in order to study its in¬ 
fluence on the cosmogenic neutrino spectrum. 

The integration procedure also has to account for the 
modification of the cosmic ray spectrum owing to inter¬ 
actions with MBR photons. This was done in two crude, 
but reasonable, ways. The first one is the introduction 
of a high energy cutoff of the spectrum as a function 
of the redshift. The second one, which is used in the 
results presented below, is to weight the cosmic ray in¬ 
jection spectrum with the interaction length Air,b on the 
infrared background radiation. The cosmic rays interact¬ 
ing in the IBR used in the integration of the yields are 
For AiRs/Atot) where Atot is the interaction length in 
the total IRB and MBR fields. The fraction of the cos¬ 
mic ray flux used in the integration procedure is shown 
in Fig. ^ If one arbitrarily determines the high energy 



The second phase of the calculation requires the 
parametrization of the redshift evolution of the emis¬ 
sivity of cosmic ray sources, and the form of the cos¬ 
mic ray injection spectrum. We assume a cosmic ray 
injection spectrum of the power-law form dN/dEp = 
exp{-Ep/E^s,yi) with E^ax = 10®i ® eV. 

We consider here two empirically based models for the 
evolution of UHECR power with redshift, viz., (I) one 
based on the redshift evolution of the star formation rate 
that was used in the calculation of the infrared back¬ 
ground in SMS05, and (2) the other based on the red- 


FIG. 4: Fraction of the total cosmic ray flux used in the in¬ 
tegration of the neutrino yields from interactions in the IRB. 
The different lines correspond to fractions at different red- 
shifts as indicated by the numbers in the plot. 


cutoff of the cosmic ray energy spectrum as the energy 
at which only 10 per cent of the cosmic rays interact in 
the IRB, it would be 7x10^° eV at 2 = 0 compared to 
2.5x101° and 4x10^® eV at 2 = 1 and 5. 

Since the yields include the dt/dz factor the integration 
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becomes very simple, viz.: 

dNl/AE^= [ dzx 

Jo 


dNp 

dEr, 


£{z)Y^[{l + z)E,;Ep,z]dEp , (3) 


where the index i indicates the neutrino flavor. 


IV. RESULTS 

The results of the integration are shown in Fig.0 The 
top panel of the figures compares the fluxes of cosmo¬ 
gonic Vp -I- Pp neutrinos generated by interactions with 
MBR photons (histogram) with those generated by in¬ 
teractions with IRB photons (squares), assuming a 7 = 1 
UHECR injection spectrum with the fast evolution of the 
emissivity of the cosmic ray sources. All panels of Fig.0 
are calculated with the same cosmological evolution. Us¬ 
ing the baseline evolution model will give neutrino fluxes 
which are about 25-30% lower. 

The peak flux of the IRB-generated neutrinos is not 
much lower than that of the MBR-generated ones, i.e., 
1.7x10“^^ compared to 2.2x10”^^ cm“^s“^sr“^. The 
peak is, however, shifted to lower by about a factor 
of 3. The main reason for that shift is the contribution 
of protons of energy below 3x10^® eV to the neutrino 
production. The IRB-generated neutrino flux is also de¬ 
pleted at energies above 10^® eV. This is because protons 
of energy above 5x10^® eV very rarely interact with IRB 
photons before they lose their energy in MBR interac¬ 
tions. At energies below the peak the IRB-generated 
neutrino flux is somewhat flatter than the MBR one, al¬ 
though the statistical uncertainty of the calculation does 
not allow us to make a quantitative statement regarding 
this. 

The middle panel of the figure shows the IRB- 
generated fluxes of Ve’s and Pg’s assuming an E~^ in¬ 
jection spectrum. The electron neutrino flux peaks at 
the same energy as the muon neutrino one. The Pe flux, 
which is due mostly to neutron decay neutrinos, is shifted 
and widened at its lower energy end. The dip between 
the iZg and Pg peaks is not as deep as it is in the MBR 
neutrino case. The reason for that is that the iZg peak is 
somewhat wider at energies below the peak. The bottom 
panel of Fig. I^shows the fluxes of IRB-generated -|- Pp 
assuming a steeper injection spectrum 7 = 1.5 and m = 
3.1. There are two main differences from the 7=1 case. 
The peak of the IRB-generated neutrino flux is higher by 
almost a factor of 3 (6.5x10“^^ in the same units) than 
for the MBR-generated neutrinos and this peak and is 
shifted down in energy by a factor of ~ 3 to ~ eV. 

This is caused by the higher flux of protons of energy 
below the high energy cutoff. In the case of the MBR- 
generated neutrinos, the general effect is not as strong 
but is reversed; the steeper proton spectrum results in a 
lower flux. 
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FIG. 5: Top panel: i/p (solid) and Pp (dashed) spectra gen¬ 
erated by interactions with IRB photons. Their sum (open 
squares) is compared to those generated by interactions with 
MBR (MBR) photons (full squares) for 7=1 assuming fast 
evolution of the cosmic ray source emissivity. Middle panel: 
Ve (solid) and Pg (dash) spectra for injection spectra as in the 
top panel. Their sum is shown with open squares. Bottom 
panel: i^p -|- Pp spectra for injection spectrum with 7 = 1 . 5 . 


Because of the very strong dependence of the flux of 
cosmogenic neutrinos on the cosmological evolution of 
the cosmic ray sources [ 2 ^, we investigated this depen¬ 
dence further. Fig. 0 compares the baseline and fast cos¬ 
mological evolutions of Ref. ^ to these of Refs. [^1^ . 
All evolution models are normalized to 1 at present, i.e. 
for z = 0 . 
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FIG. 6 : Four different models for the cosmological evolu¬ 
tion of the cosmic ray sources (see text for description and 
references). 


FIG. 7: Cosmogenic nentrino fluxes for four different models 
of the cosmological evolution of the cosmic ray sources - see 
text. 


The evolution taken from Ref. |3, that was used for 
calculations of the cosmogenic neutrino flux from interac¬ 
tions in the MBR 0 , gives a UHECR emissivity which 
is about 60% higher at redshift of 2 than the average of 
the models of SMS05. The cosmological evolution of the 
flat spectrum radio galaxies |2l| has an intermediate red¬ 
shift evolution; it is faster than m = 3 and slower than 
TO = 4 below z = 1 and peaks at about z = 2. It is also 
distinguished by its rapid decrease in emissivity at z > 3. 

Figure [7| compares the cosmogenic neutrino fluxes of 
generated by the baseline and fast models of 
SMS05 with the to = 3 model used in Ref. that 

of Ref. . The difference in the calculated fluxes is ac¬ 
tually quite small, compared to all other uncertainties of 
the calculation. The fast and the baseline models bracket 
from above and from below the fluxes of cosmogenic neu¬ 
trinos from interactions in the IRB, while the other two 
models fall between these two. The main reason for the 
small differences is the dz/dt matrix element that de¬ 
creases the contribution of higher redshifts because in 
the cosmological integration the emissivity is multiplied 
by the smaller time intervals involved at higher redshifts. 

In the case where UHECR luminosity evolution is as¬ 
sumed to be proportional to the redshift distribution of 
flat spectrum radio galaxies as given in Ref. , the neu¬ 
trino spectra peak at a somewhat higher energy because 
of the relatively small UHECR emissivity at the higher 
redshifts. 

In Fig.|Slwe present the total fluxes of cosmogenic neu¬ 
trinos from interactions in the MBR and IRB for injection 
spectra with 7 = 1.0 and 1.5 and for fast cosmological 
evolution of the cosmic ray sources as in Ref. 0 - 

For a relatively flat (7 = 1) injection spectrum (empty 
squares in Fig. 0 ) interactions with IRB photons gener¬ 
ate almost as many cosmogenic neutrinos as interactions 
on MBR. The peak of the total neutrino energy spec¬ 
trum from interactions in the MBR and in IRB is shifted 
to lower energy by a small amount (from 3x10^^ eV to 



FIG. 8 : Total muon + fluxes for 7 = 1.0 (empty squares) 
and 1.5 (fnll squares) calculated with fast cosmic ray source 
evolution. The shaded area represents the W&B 1^ upper 
bound on astrophysical source neutrinos. 

about 2x 10^^ eV). The distribution extends to lower neu¬ 
trino energies by more than half a decade. 

For steeper spectra (7 = 1.5) the contribution of IRB- 
generated neutrinos is more significant and the result¬ 
ing flux is almost an order of magnitude larger. The 
magnitude of the flux at the peak of the spectrum is 
~ 10 “^® cm“^s“^sr“^ and is higher than that of the 
MBR-generated neutrinos by a factor of ^ 3. 

In both cases there is no increase of the neutrino flux 
at energies exceeding 10^® eV. The influence of increased 
cosmological evolution of the cosmic ray sources is practi¬ 
cally the same as in the case of MBR-generated neutrinos 
alone. In this context, we note that the cosmological evo¬ 
lution of target IR photon density is slower that of the 
MBR. 

For comparison, the shaded area in Fig. 0 shows the 
upper bound on the astrophysical neutrino spectra given 
in Ref. |^. The lower edge is the bound in absence of 
cosmic ray source cosmological evolution, and the upper 
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TABLE I: Rates per km® water per year of showers above 
different energy generated by different types of neutrino in¬ 
teractions for cosmic ray power density Vo at z=0 of 1.4x 10®® 
W Mpc“® and fast cosmological evolution for homogeneously 
distributed cosmic ray sources (see text). 


logio Esh 

7" 

=1 

7= 

1.5 

(GeV) > 

MBR 

IRB 

MBR 

IRB 

6 

0.092 

0.021 

0.078 

0.085 

7 

0.088 

0.019 

0.072 

0.072 

8 

0.079 

0.010 

0.063 

0.030 

9 

0.044 

0.001 

0.027 

0.001 

10 

0.008 

0.000 

0.003 

0.000 


edge is for (1 + z)® evolution. 

The inclusion of the proton-IRB interactions somewhat 
reverses the trend of the injection spectrum dependence 
of the cosmogenic neutrino flux. Without including such 
interactions, steeper injection spectra lead to smaller cos¬ 
mogenic neutrino fluxes; with the inclusion of the contri¬ 
bution to the neutrino flux from interactions of IRB pho¬ 
tons with relatively lower energy protons, steeper cosmic 
ray spectra generate higher neutrino fluxes. The reason 
is that we normalize the cosmic ray injection spectrum at 
10®® eV, which is now in the middle of the energy range 
of the interacting cosmic rays. One can see in Fig.^jthe 
dominance of the interactions in the MBR of cosmic rays 
of energy above 10®® at all redshifts higher than 1. For 
steeper cosmic ray injection spectra the number of such 
particles is decreased while that of cosmic rays below 10 ®® 
eV, that interact in the IRB, is significantly increased. 

Because of the lower average energy of the IRB- 
generated neutrinos, the spectra are shifted and their de¬ 
tectability is lower than that of MBR generated photons. 
This is because the neutrino-nucleon cross section rises 
monotonically with energy. Table Q] shows the shower 
rates of and 9e CC (charged current) interactions per 
km®yr of water detector for cosmogenic neutrinos gener¬ 
ated by interactions with MBR and IRB photons. These 
rates are the products of the neutrino cross section times 
the neutrino flux integrated above Egh- We assume that 
the total neutrino energy is transfered to the shower ini¬ 
tiated by its CC interaction. The calculation of event 
rates of CC and NC interactions of muon and tau neu¬ 
trinos are much more difficult and require Monte Carlo 
models of particular experiments. 

The second column of Table U corresponds to the num¬ 
bers given in a similar table in Ref. |l2l . The numbers 
can not be directly compared because of the different cos¬ 
mologies (LIm = 1 in Ref. and LIm = 0.3 here) and 
cosmological evolutions of the cosmic ray sources used. 
The = 0.3 cosmology increases the neutrino rates by 
about 70%. 

In the 7 = 1.0 injection case the rates from MBR neu¬ 
trinos are higher by factors above 4 at all shower thresh¬ 
olds, while in the 7 = 1.5 case the IRB rates are higher 
or similar for shower thresholds below 10® GeV. Above 
that energy MBR neutrinos generate higher rate. Note 



FIG. 9: Total muon i/i_i + Vp, fluxes for 7 = 1.7 and evolution 
from Ref. source evolution (full squares) and no evolution 
(empty squares). 


that IRB neutrinos do not contribute at all to the shower 
rates above 10® GeV. 

The total shower rate for Egh > 10® GeV is higher than 
the MBR only rate by about 20% in the 7 = 1.0 case, while 
in the 7 = 1.5 case it more than doubles the detection 
rate. 

Most of the contemporary fits of the injection spectrum 
of the highest energy cosmic rays confirm the conclusion 
of Ref. that it is steeper than The spectrum 

derived by these authors is E“® ®’ with a signihcant flat¬ 
tening at about 10®® eV, which could be explained with 
different effects, see e.g. Ref. The shape of the spec¬ 
trum may be accounted for in this model as a result of the 
P 7 ^ e+e“ process as discussed in Ref. [^. This 
pair production process creates a dip at about 10®® eV 
and a slight excess at the transition from pair production 
to purely adiabatic proton energy loss at about 10®® eV. 
This ht does not require a strong cosmological evolution 
of the cosmic ray sources, but can accommodate a mild 
one (X (1 -I- z)"* with m < 3 113 ■ In the case of flatter 
injection spectrum the pair production dip is well ht also 
by TO = 4. Figure El shows the spectra of the cosmogenic 
neutrinos from interactions with MBR and IRB photons 
assuming a steep injection spectrum with 7 = 1.7 and (1) 
no evolution with (to = 0 ) and ( 2 ) evolution according to 
Ref. [IJ. The difference between the two neutrino spec¬ 
tra is signihcant; the peak values are 10“®®"(1.5xl0“®®) 
cm“®s“®sr“® for without and with evolution. The ad¬ 
dition of the IRB component brings these spectra into 
the range of detectability, especially in the case of mild 
cosmological evolution. 


V. DISCUSSION AND CONCLUSIONS 

The evolution models of SMS05 do not give the highest 
IRB-generated neutrino hux. We compared the IRB pho¬ 
ton density in this model with the models of Refs. p^l34l| . 








7 


Both of these models have higher IRB density in the rel¬ 
evant energy range between 3x10“^ and 1 eV. The total 
IRB densities in the this range are 1.27 and 1.12 
compared with the density of 1.03 used here. In addi¬ 
tion, Ref. shows somewhat faster cosmological evo¬ 
lution. The use of any of these models would have in¬ 
creased somewhat the calculated flux of cosmogenic neu¬ 
trinos. The uncertainty in the IRB flux is of the order of 
30% , while the biggest uncertainty in this calculation 

is in the UHECR flux and its cosmological evolution. 

The IRB contribution to the total cosmogenic neutri¬ 
nos flux can be slightly increased as protons of energy 
below 10^® eV can interact with IRB photons and gener¬ 
ate lower energy neutrinos. If such interactions were in¬ 
cluded the IRB spectrum would be wider than the MBR 
one, especially at energies below 10^® eV. 

It is difficult to compare our results with those of 
Refs. [I1II3 because of the different astrophysical input 
in these calculations. Qualitatively the results of these 
calculations are similar to ours and certainly agree within 
a factor of 2 . 

In conclusion, we calculated the flux of cosmogenic 
neutrinos from interactions of UHECR protons with IRB 


photons using the recent calculations of IR photon spec¬ 
tra densities as a function of redshift by SMS05. Our 
calculations show that UHECR interactions with IRB 
photons produce a signihcant flux of cosmogenic neutri¬ 
nos, one which is comparable to interactions with MBR 
photons. This is especially true in the case of assumed 
steep injection spectra of the ultrahigh energy cosmic 
rays. The total neutrino event rates at energies above 
I PeV increase by more than a factor of 2 in the case 
of injection spectra with 7 = 1.5. Because of the much 
lower mean free path of protons above I0^° eV in the 
MBR interactions with IRB photons do not increase the 
higher energy end of the cosmogenic neutrino spectrum. 
The total cosmogenic fluxes, however, are still not de¬ 
tectable with conventional neutrino telesco pes such as 
IceCube or the European km® telescope J 3b | . A reli¬ 
able detection is only expected from radio and acous¬ 
tic neutrino detectors. 
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